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EQUATIONS  OF  MOTION  FOR  A  TOWED  BODY 
MOVING  IN  A  VERTICAL  PLANE 


INTRODUCTION 

The  designer  of  a  ship- towed  system  must  be  able  to  predict 
the  motions  of  the  towed  body.  In  designing  a  body  that  should 
have  minimal  response  to  external  excitations,  he  should  have 
such  information  available  in  order  to  avoid  a  resonance  in  any 
of  the  body  motion  modes. 

This  report  presents  solutions  of  the  cross-coupled  equations 
of  motion  for  a  towed  body  moving  with  heave,  surge,  and  pitch 
motions  and  for  various  ship-input  conditions.  By  the  use  of  body 
and  cable  constants,  one  can  optimize  a  design  for  whichever 
criterion  he  must  satisfy,  be  it  heave,  surge,  or  pitch,  or  any 
combination  of  the  three  motions.  The  analysis  is  limited  to 
motions  in  the  vertical  plane,  since  additional  information  that 
could  be  derived  from  solving  equations  of  motion  for  a  body  with 
six  degrees  of  freedom  did  not  appear  to  be  significant  at  the  time. 
In  general,  towed  bodies  have  symmetry  about  a  plane  defined  by 
the  vertical  and  the  fore-and-aft  axes. 

The  primary  force  input  to  a  cable-towed  system  is  taken  to 
be  the  vertical  force  at  the  point  of  attachment  of  the  cable  on  the 
ship.  This  limitation  on  the  analysis  is  reasonable,  because 
forces  in  this  plane  will  not  couple  into  sway,  roll,  or  yaw  for  a 
symmetrical  body.  The  solution  may  be  extended,  however,  to 
apply  to  a  body  moving  with  all  six  degrees  of  freedom. 

Lum1  solves  similar  equations  by  using  the  Laplace  transform 
technique  for  the  open-loop  pitch  response  of  a  towed  body;  how¬ 
ever,  the  results  are  somewhat  difficult  to  use  to  find  towed-body 
response  for  specific  input  conditions. 


*S.  M.  Y.  Lum,  A  Theoretical  Investigation  of  the  Body  Parameters  Affecting  the  Open- 
Loop  Pitch  Response  of  a  Submerged  Towed  Body,  David  Taylor  Model  Basin  Report  No.  1369, 
February  1960  (UNCLASSIFIED). 


ANALYSIS 


A  typical  towed  system  is  shown  in  Fig.  I.  The  dynamic 
model  of  the  system  is  presented  in  Fig.  2. 


Fig.  1  -  A  Tewed  System 

The  dimensions,  masses,  spring  constants,  and  damping  fac¬ 
tors  for  Fig.  2  are  defined  as  follows: 

Cfay  -  damping  factor  of  the  body  alone  for  motion  in  the  y-direction 

Chz  -  damping  factor  of  the  body  alone  for  motion  in  the  z-direction 

Cjjr  -  rotational  damping  factor  of  the  body  alone  for  motion  about  the  x*axis 

C{  -  damping  factor  of  the  tail  alone  for  motion  in  the  z-dircction 

D  -  vertical  distance  from  center  of  mass  to  tail  center  of  damping 

G  -  vertical  distance  from  the  center  of  mass  vo  die  tail's  effective  center 

of  resilience  and  damping. 

h  -  vertical  distance  from  the  center  of  mass  to  the  towpoint 

j  -  vertical  distance  from  the  center  of  mass  to  the  body’s  center  of  damping 

Kb  -  spring  constant  due  to  change  of  body  lift  with  angle  of  attack 
Kcy  -  spring  constant  of  deflecter  towcable  in  horizontal  (y-direction) 
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Fig.  2  -  Dynamic  Model  of  die  Towed  System 
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-  spring  constant  of  deflected  towcable  in  vertical  (z-direction) 

-  spring  constant  due  to  change  of  tail  lift  with  angle  of  attack 

-  sj.  ing  constant  due  to  change  of  tail  drag  with  angle  of  attack 

-  rotational  spring  constant  due  to  pendulous  effect  of  the  water  weight  of 
the  body  below  the  towpoint 

•  horizontal  distance  from  front  of  body  to  the  center  of  body  lift 

-  horizontal  distance  from  front  of  body  to  the  center  of  body  damping 

-  horizontal  distance  from  front  of  body  to  the  center  of  mass 

-  horizontal  distance  from  front  of  body  to  the  effective  center  of  the  tail 

-  horizontal  distance  from  front  of  body  to  the  towpoint 

-  mass  of  the  body  and  tail 
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mass  of  the  towcable 

mass  of  vater  that  floods  body 

hydrodynamic  mass  of  the  body  ic  the  y-directioc 

hydrodynamic  mass  of  the  body  in  the  z-direction 

hydrodynamic  mass  of  the  cable  in  the  y-direction 

hydrodynamic  mass  of  the  cable  in  the  z-direction 

hydrodynamic  mass  of  the  tail  in  the  z-direction 

vertical  distance  from  towpoint  to  center  of  gravity  of  the  dry  body 

vertical  displacement  of  the  towpoint  on  the  ship 

weight  of  the  body  in  sea  water 

horizontal  distance  from  front  of  body  to  center  of  gravity  of  the  dry  body 
horizontal  distance  from  front  of  body  to  center  of  body  hydrodynamic  mass 
horizontal  distance  from  front  of  body  to  center  of  body  flood  water  mass 


The  center  of  mass  of  the  body  can  be  located  by  taking 
moments  about  some  reference  point.  For  its  location  in  the 
y-direction 


'CM 


mb +  yj  mhbz +  y3  ®fw +  Lc  mht +  ltp  (1/?  mc + 1/3  mhcz> 

“b  +  “hbz  +  rafw  +  “hr  +  1/3  mc  +  */3  mhcz 


(1) 


The  mass  and  hydrodynamic  mass  terms  of  the  cable  are  included 
in  Eq.  (1),  following  the  usual  convention  for  dynamic  problems 
to  include  one  third  of  the  mass  of  the  spring.  A  more  thorough 
analysis  might  be  patterned  after  the  work  of  Dr.  L.  F.  Whicker.2 

The  rotational  spring  constant  Kw  is  equal  to  wbn  if 
angular  displacements  are  assumed  to  be  small.  The  spring 
constants  Kb  and  Kt  are  actually  rotational  spring  constants 
but  are  shown  as  linear  springs  in  the  z-direction  for  ease  of 
computation.  The  constants  can  be  computed  as  follows: 


2L.  F,  Whicker,  "The  Oscillatory  Motion  of  Cable-Towed  Bodies,"  Doctoral  Thesis, 

University  of  California,  Berkeley,  California,  May  1957, 

<r 
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Kb  =  Vi2  AbV2(AClb), 


(2) 


Kt  =  />/2AfV2(AClt),  (3) 

and 

KDt  =  p/2A,V2(ACDt),  (4) 

where 

p  =  fluid  density, 

Ab  =  horizontal  projected  area  of  the  body, 

A(  =  horizontal  projected  area  of  the  tail, 

AC|b  --  change  in  lift  coefficient  of  the  body  with  angle  of  attack, 

ACm  =  change  in  lift  coefficient  of  the  body  with  angle  of  attack,  and 

A CQt  =  change  of  drag  coefficient  of  the  tail  with  angle  of  attack. 


A  criterion  for  stability  is  that  the  lift  moment  of  the  tail  must 
be  greater  than  the  lift  moment  of  the  body;  it  can  be  expressed 
mathematically  by 


Kt(Lt-LCM)  >  %  U-cm- l‘bi 


Using  the  equations  of  motion  for  a  resiliently  supported  body3 
with  six  degrees  of  freedom  as  presented  by  Harris  and  Crede,4 


^This  form  of  equation  of  motion  is  used  because  of  its  linear  form.  Viscous  forces  are 
represented  as  damping  constants  and  inertial  history  terms  are  assumed  to  be  negligible, 

4C.  M.  Harris  and  C.  E.  Crede,  Shock  and  Vibration  Handbook,  vol.  1,  McGraw-Hill 
Book  Company,  iuc.  ,  New  York,  1961, 
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but  restricting  the  motion  to  a  plane  with  inputs  occurring  in  that 
plane  only,  we  can  write 
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zzAy“Kyz  Az)(Zc-W) 
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Al  *  AJ  -  A,  A.)  “  +  2 (C"  Ar'  _tVy  A«>  *' 

A,'  -  c,y  Ai>  *-*>♦£  (cyy  A»!  +  C=  V  *  2Cy»  Ay  A  V  « • 


(5) 


FV  *  "y  +£  Kyy  XI  K>z  <I«-*>+£  (Ky*  Ay  'Kyy  A*>  “ 
*  £Cyy  ^  +  £  S*  *-**♦£««>.  -  Cyy  AP  » • 


(6) 


and 


Fz  “  mz  zc  +  yi  ^yz  yC  +  ^  1  Kzz  <zc*  *>  +  <Kzz  Ay  Kyz  Az)  “ 

*  £  Cy*  1  £  c,z  -  »>  +  S  <c“  Ay  -  cy«  Az»  “• 


(7) 


where 

A  «  coordinate  distance  from  the  body’s  center  of  mass  to  the  elastic  center 
of  the  resilient  element; 

A'  *=  coordinate  distance  from  the  body’s  center  of  mass  to  the  effective  center 
of  the  damping  element; 

C  e  damping  factors; 

F  =  dynamic  input  force  applied  at  the  body’s  center  of  mass; 

I  «  mass  movement  of  inertia; 

K  •=  spring  constant; 
m  =  body  mass; 

M  *  dynamic  input  moment  applied  at  the  body’s  center  of  mass; 

W  =  vertical  displacement  of  the  ground  point  of  all  vertical  springs; 
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V,  y,  y  =  displacement,  velocity,  and  acceleration  in  horizontal  direction  of  the 
body's  center  of  mass  relative  to  earth  coordinates; 

z,  z,  z  -  displacement,  velocity,  and  acceleration  in  the  vertical  direction 
of  the  body’s  center  of  mass  relative  to  earth  coordinates; 
a,  a,  a  =>  angular  displacement,  velocity,  and  acceleration  in  a  vertical  plane 
(defined  by  y  and  z)  of  the  body  relative  to  earth  coordinates; 


and  where 

subscript  c  =• 
subscript  x  = 
subscript  y  = 
subscript  z  = 
subscript  a  = 


reference  to  center  of  mass; 
reference  to  x-direction; 
reference  to  y-direction; 
reference  to  z-direction;  and 
reference  to  a -axis  of  rotation  . 


The  computation  of  the  various  coordinates  distances  A  and 
A'  is  as  follows: 

for  the  spring  representing  the  cable, 

Ay  =  lcm  “  ltp 
Az  =  h  ; 


for  the  spring  due  to  tail  lift, 


Ay  “  lcm  "  Lc 


Az  =  G  (negative  if  effective  center  of  tail  lift  is  below  the  center  of  mass); 


for  the  spring  due  to  body  lift, 

Ay  =  LCM  ~  Ht 

A?  =  0  (it  is  assumed  that  the  effective  center  of  body  lift  lies  on 
the  horizontal  plane  of  the  center  of  mass); 


•mi 
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for  the  center  of  damping  of  the  body. 


“lcd 


K  ~  J  (eegacire  if  tbe  center  of  dumping  lies  below 


*e  center  of  mu}; 


for  the  center  of  damping  of  the  tail. 


Ay  *  Lcw  ~Lt 


A z  *  G  (negative  if  die  center  of  tail  dating  is  below  the  center  of  amss). 


We  may  now  consider  the  actual  towed  system  and  use  all 
masses,  elastic  elements,  damping  elements,  and  input  forces 

/At°r  app!y  and  (7).  P.ewriting  Eqs.  (5), 

),  and  (7)  with  the  above  conditions  and  coordinate  distances 
and  assuming  that 


Kyz  *  KZy  =  0  for  tfre  cable,  and 
Cyz  *  C*y  =  0  for  the  body. 


we  have 


M*  =  ^b+/hb>°+(“Kchh>yc+tKcv(LCM-Llp)](zc-W) 

+  tKch  h  +  Kcv  (Lcm  ~  Ltp>2  +  K,  <L,  "  hCM)  -e  Kb  (Lcm  -  Lbl)  +  Kw  +  Kdt  D]  a 
-^bh  i)  yc  +  ICbv  (Lcm  ~  Lc  d)  -  c,  (Lt  -  lcm)  j  ic 
+  tCbh  i 2  +  C,  fL,  -  Lcm)2  +  Cby  (Lcm  -  Lco)2  +  Cbr  I  a. 


(9) 


Fy  *  "«y  >c  4  Kci,  >c  "  (KcS,  fc  "  Kd«>  °  4  ^  4  lCb/Lb  ~  Ct*  >j 


and 


F,  =  f.  *  K„  (zc  -  *)  *  |I„  (Lcu  -  Ltp)  *  Kt  -  K,  I  a 

~  K  4  ^“Ci  *Lt  ”  LOi*  4  ^LCM  “LCD^  4  a  > 


(10) 


where 
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“«y 


ra 


tz 


•ass  soacfli  of  inertia  in  the  pirdi  plane  about  the  center  of  aass; 

hydrodynamic  mass  moment  of  inertia  in  tbe  pitch  plane  about  tbe  center  of  mass; 

forcing  function  acting  on  tbe  center  of  aass  in  tbe  y-directicw; 

forcing  function  acting  on  tbe  center  of  mass  in  tbe  z -direct!  on; 

forcing  function  acting  about  die  center  of  mass  in  tbe  pitch  plane; 

total  oscillating  aass  in  the  y-direction  (mfy  -  mft  +  mhby  +  ®f  w 
+  1/3  mc  +  1/3  mhcy);  and 

total  oscillating  aass  in  the  z-direction  (mtz  -  ab  +•  *j,5z  +  m{  w 
*  mht  4  ’/3  “c  4  1/3«hc*>- 


The  spring  constant  due  to  body  lift  is  negative  in  Eq.  (8) 
because  its  effect  is  to  upset  the  body  if  it  is  displaced.  Other 
assumptions  applied  to  Eqs.  (8),  (9),  and  (10)  are  that  cable 
damping  is  small  compared  with  body  damping  and  that  the 
center  of  tail  damping  lies  along  the  vertical  height  of  the  center 
of  mass.  It  is  also  assumed  that  the  flood  water  mass  center, 
the  hydrodynamic  mass  center,  and  the  center  of  body  heave 
damping  lie  in  the  same  horizontal  plane. 
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SOLUTION  OF  EQUATIONS 


In  order  to  change  the  equations  of  motion  (Eqs.  (8),  (9),  and 
(10)  )  into  dimensionless  form,  the  following  steps  are  taken: 

(1)  All  a  terms  are  multiplied  by  Lfc/Lb,  where  Lb  is 
the  body  length. 

(2)  The  two  force  equations  are  divided  by  wb  ,  the  body 
weight  in  water. 

(3)  The  moment  equation  is  divided  by  wb  Lj, . 

Therefore,  Eqs.  (8),  (9),  and  (10)  become 


!« +Kcr*Lcn  ~ltp)W  Jb  +  Jhfc  _  M  P^cs-Wp) 

Lb*b  .  Lb  wb .  b  Lb  *b  fc  rTiT  *c 


TKch  +  Kcv*LCM  ~  LTP^  +  Kc<Lt  “  + 

L  ^ 


Kb^LCU  **Lb|) 


^Kw^KdcDj 


f  cbh  J]  .  fo*  <Lc*  ~Lco>  +Ct<Lc«  -Lt)‘  . 

'Mn — ^ — n* 


(aLb) 

(I) 


Cbh  r  *  Ct  ( LCU  -  Lt>2  +  Cbr  <LC«  -  LCD>2  +  Cbr 


Lb  Wb 


(Lb  a). 


(II) 
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{Ill) 


For  Eqs.  (I),  (II),  and  (III),  let 


Lb  *  =  q, 
Lfc  a  =  4  , 


and 

* 


Lba  =  q. 

In  the  equations  that  follow,  the  terms  I,  II,  and  III  indicate 
Eqs.  (I),  (II),  and  (III),  respectively,  and  the  subscript 
following  these  terms  represents  the  variable  that  multiplies 
the  term. 


I 


Kk  (LCB  -  Ltp)!  *  K,  (L,  -  Lcm)  +  Kb  fLcu  -  Lbl)  *K,tKJtD‘ 

Lb  *b 


fCbh  j  Cbr 

II  •  *=  !  - -  +  - 

q  Wb  Lb  tv  L2 

L  b 


* 

* 


i 


Z  =  III“i+IIIzi  +  IIIqq  flll^z+III^q;  (I HA) 

or  in  matrix  notation. 


or 

MQ(t)  +  C()(t)  rKQ(t)  =  F(t).  (11) 


The  solution  to  Eq.  (11)  involves  a  complementary  and  a 
particular  integral.  The  complementary  solution,  which  is 
obtained  by  setting  the  forcing  function  F{t)  of  Eq.  (11) 
equal  to  zero,  is 

Qc  “  <-'(C/zm)t  (C,  cos  c,nd{  +  C2  sin  ojndt), 


V- 


where 

r  =  time, 

Cj  an d  C2  =  arbitrary  constants,  and 

o>n(j  -  damped  natural  frequency  of  tbe  system. 

It  can  be  seen  that  the  complementary  solution  is  a  transient  one 
that  diminishes  with  time.  Therefore,  the  steady-state  condition 
is  described  by  the  particular  solution. 

To  effect  the  particular  solution  we  let 

F  (t)  =  F0  cos  Git , 


and  assume  that 

0  =  Q2  cos  rut  +  Q j  sin  (ot, 

where 

a  =  forcing  frequency. 


Taking  derivatives  and  substituting  in  Eq.  (11)  we  obtain 


0  =-a  Q2  sin  tit  +  a  Qj  cos  at  , 


and 


0  *=  -  a2  Q,  cos  at  -  aA  Qj  sin  at  J 


therefore,  we  have 


-gj2  M  (Q2  cos  gh  +  Q  j  sin  gh)  +  gj  C(Qj  cos  gh  -  Q2  sin  Git) 


+  K(Q2  cos  at  H  Qj  sin  Git)  =  F0  cos  gh 


(12) 
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Equation  (12)  must  satisfy  the  condition  that  the  sum  of  the  forces 
in  the  vertical  and  horizontal  directions  must  equal  the  inertia 
forces  in  those  directions.  Hence,  we  have 


-  oj2  MQ2  +  ojCQ1  +  KQ2  =  F, 


-  6i2  MQt  -  &>CQ2  +  KQt  =  0.  (; 

Because  M,  C,  and  K  are  matrices,  multiplication  of  these 
terms  is  not  commutative. 

Solving  Eq.  (14)  for  Q2  and  substituting  into  Eq,  (13), 
we  obtain 


Q2  =  JLc*1  l  K  -  a)2  M  ]  Qj  , 


j  i  C* 1  (K  -  oj2  M)  ic'!(K-w2M)+  1  JO,  =  -ic‘J  F. 
L  to  (o  J  a) 


Solving  for  Q},  we  have 


Qj  *  ric*l(K-w2M)2  +  ll’1  ic^F, 

U  J  (O 


in  which 


"y  0  -H<i 


C  -  0  III^  III^ 


- 1  •  1  *  I  • 

y  z  q 


and 
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ic*  (K  -  M)  .1  jc^fl Ir  -  J  IIf)  .  IyCy,  Cy,ail,  -  III;)  +  C„ I, 
-C,«  *c„  «  C„0,  -  J  I Jj)  Cly  (IIy  -  c*  II-) 

“IyC2q  ^q^  +  Czx  (IIIZ-6>2lIIi)-Czynq  +  Cxx!II<i 

+  Czq«q  ~  "2  Iq>  Cqy(Hy  -  a>2  II-)  -  CqqIy  Cqx  (III2  -  o>2  III-) 

+  c«  '•  •c,",'c»111,*  c«  c,  -  ■»*  I,)]  • 
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Now  forming,  we  have 
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rqq 
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rqy 

rqz 
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lid#*  *  - - -  -  *  ~ 


(*yy  4  ryz  rzy 

ryy  *y*  +(u  [y* 

*ryy  ryq  +  ry*r*q 

+  r<lqrqy)(1/<a2)  +  1 

+  rqqryq>(l/w2) 

^rzy  ryy  4  rzz  rzy 

<rzyryz+r« 

'ryq‘*y  +  fzqr*z 

+  rzqrqy)(I/6,2) 

+  rqzrzq^1/6,2)+  1 

+  fzq  rqq>  <l/t,2) 

^rqy  ryy  +  rq*  rzy 

^fyz  fqy  4  fzz  fqz 

+  rqqrqz)(l/cj2) 

+  rqzrzq)  d/<02)  +  1 

Let 


Forming,  we  have 


[(ft»C*!  (K.  -  Cl)2  M))2  +  l  I'1  = 
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where 


J„  '  VB- 

^*y 

^qy  *  b*ybq*“bKbqj/®» 

*y*  =  hyql,q*“by*hq/B  * 

J«  “  hyy h<q“byqhq/B* 

^q*  by*bqy  ~byy  bqz^B  * 

Jy,  "  hyxi,*q-hzzhy«/B* 

JZq  “  hyqb*y  “hyy  ^zq^®* 

^qq  *  Sy8M~8y*8*/8, 

where 

B  “  hyy(h«hqq"hzqhqz)+hy*(hzqhqy”!,zyh,^+hyq(hzyhql:-»'xs 
Forming,  we  have 


Ic*1  F 

m  J_ 

c 

yy 

c 

s* 

c 

Cyq 

c 

CD 

o 

*y 

zz 

*q 

-Cqy 

C 

q* 

c« 

and 


CyyV+CyiZ+Cy,p' 

r  - 

E 

clyV  +  C,£ZtCl<1P 

1 

TZ 

Cl) 

F 

.c1> * + c,<  ~z>  c„, p. 

G 
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,i  m  <\Jj 

-?*"*** 


Z3G&WZ5ZZX 


Finally,  v/e  obtain 


and 


-  - 

- 

*2 

9n*t 

ryz  z  *  + 

1 

Cl 

*  r«Z!  "  'z,*! 

q* 

*  '«**«  * 

y  »  y ,  sin  ioz  *  y2  cos  ot, 

z  »  2j  sin  £j r  +  i,  cos  ot,  and 

q  =  qj  sin  <jt  4  q,  cos  wl 


COMMENTS 

A  simple  computer  program  has  been  written  to  effect  these 
complex  solutions.  This  program  with  necessary  inputs  is  shown 
in  Appendix  A. 

in  working  with  these  equations,  one  finds  that  the  most  diffi¬ 
cult  task  is  determining  the  cable  spring  constants,  hydrodynamic 
masses,  and  damping  factors  since  there  is  very  little  information 
in  the  literature  dealing  with  the  computation  of  these  factors  for 
towed  bodies.  Appendix  B  discusses  the  computations. 

The  solutions  have  been  used  to  study  the  motions  of  the 
AN/SQA-10  towed  body  for  various  ship  speeds  and  sea  conditions. 
The  results  agree  favorably  with  data  taken  at  Sea  States  2  and 
5  with  a  full-size  AN/SQA-10  towed  body. 


r, 

22 


5 


trm 


SUMMARY 


The  differential  equations  of  motion  have  been  written  for  the 
case  of  a  cable-supported  towed  body.  The  solution  for  a  sinusoidal 
forcing  function  input  is  shown  in  £qs.  (15)  and  (16).  Computa¬ 
tion  of  the  inertial  terms,  spring  constants,  and  damping  factors 
required  for  Eqs.  (12),  (13),  (14),  (15),  and  (16)  enables  one  to 
solve  for  the  magnitude  and  phase  angle  of  motions  in  surge,  heave, 
and  pitch  resulting  from  vertical  displacements  of  the  cable  end  on 
the  ship.  A  trial-and-error  approach  can  be  used  to  determine  the 
natural  frequencies  of  the  system  by  substituting  various  frequen¬ 
cies  into  the  solution. 
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Appendix  A 

COMPUTER  PROGRAM  FOR  COMPUTATION 
OF  TOWED  BODY  MOTIONS  IN  SURGE, 
HEAVE,  AND  PITCH  FOR  A 
VERTICAL  SINUSOIDAL  INPUT  AT  THE  SHIP 


INPUT  TO  COMPUTE  SURGE,  HEAVE,  AND  PITCH  OF  A 

TOWED  BODY  * 

«b  -  mass  of  body  and  tail 
Wb>  •  weight  (includes  tail)  in  air 

a^  -  hydrodynamic  mass  of  body  in  fore-and-aft  direction 
nhbz  *  hydrodynamic  mass  of  body  in  vertical  direction 
t*fw  flood  water  mass 

-  hydrodynamic  mass  of  tail  in  vertical  direction 
Vb  -  body  weight  (includes  tail)  in  water 

Vj  -  distance  from  front  of  body  to  body  center  of  gravity  in  air 

y2  -  distance  from  front  of  body  to  body  crater  of  vertical  hydrodynamic  mass 

y,  -  distance  from  front  of  body  to  flood  water  mass  center  of  gravity 

Lt  -  distance  from  front  of  body  to  tail  hydrodynamic  mass  center  of  gravity 
and  tail  damping 

LTp  -  distance  from  front  of  body  to  towpoint 

Lco  -  distance  from  front  of  body  to  body  center  of  damping 

Lbg  -  distance  from  front  of  body  to  center  of  body  lift 

n  •  distance  from  towpoint  to  body  center  of  gravity  in  air 

n  •  distance  from  towpoint  to  body  center  of  mass  in  water  (positive  if  above 
center  of  mass) 

j  -  distance  from  center  of  mass  to  center  of  damping  (minus  if  center  of  damping 
below  center  of  mass) 


*  Units  must  be  kept  homogeneous. 


Kch  *  horizontal  cable  spring  constant 
Kct  -  vertical  cable  spring  constant 


Kb 

•  body  lift  spring  constant 

C,b  p/2  Ab  V' 

=  f<«> 

•  tail  lift  spring  constant 

C„p/2A,VJ 

C|t  -  ««> 

K* 

•  tail  drag  spring  coastact 

C*  p/1  A,  W 

Cdt  =  Ha) 

V 

displacement  amplitude  of  towpoint  on  ship 

cj  •  angular  frequency  of  ship’s  pitch  motion 
J  b  -  mass  moment  of  inertia  of  body  about  center  of  mass 

-  hydrodynamic  mass  moment  of  inertia  of  body  about  center  of  mass 
Cj  -  vertical  damping  factor  of  body  alone 

Cbh  -  horizontal  damping  factor  of  body  alone 
Cbr  *  rotational  damping  factor  of  body  alone 

-  vertical  damping  factor  of  tail 

■  »c  -  mass  of  towcable 

mhcz  -  vertical  hydrodynamic  mass  of  cable 
mjjcy  -  horizcntal  hydrodynamic  mass  of  cable 


PRINTOUT 

The  surge,  heave,  and  pitch  will  be  printed  out  in  the  following 
order. 

Real  part  of  surge 
Imaginary  part  of  surge 

Vector  sum  of  real  and  imaginary  parts  of  surge 
Phase  angle  between  real  and  imaginary  parts  of  surge 
Real  part  of  heave 
Imaginary  part  of  heave 

Vector  sum  of  real  and  imaginary  parts  of  heave 
Phase  angle  between  real  and  imaginary  parts  of  heave 


Real  part  of  pitch 
Imaginary  part  of  pitch 

Vector  sum  of  real  and  imaginary  parts  of  pitch 
Phase  angle  between  real  and  imaginary  parts  of  pitch. 

The  uni; of  surge  and  pitch  will  be  commensurate  with  the 
units  of  the  input  constants.  The  phase  angles  are  in  degrees.  The 
pitch  angles  are  in  radians  times  the  body  length. 


COMPUTER  PROGRAM 

r  FO  att^as  of  motto*  for  a  rn.^pr  orry  (-.4^4  pa<pp 

5  ,iOn  'VTV,onPM'',fl4*n 
«  TNP  ^CT5'.6r AP(Sn60ft0 

m-F'  sir'1  t tyoh< 20 j  .tt*ov/<70)  ,ttypi«i  ,ttyt  t4)  .tonfomai  ,ttmcp5J  . 

»TT«'HrvfS) .  tim ?o)  ,  i«M70)  .rssdo  .  ine?o)  ,100*170)  ,sf<2o>  * 

2  (7fl)  »WTU)  •  T< f' T  f4)  ,F*'T  (]?| 

OF AO  INPUT  TAPF  3,100.1  JPR.Mnm 

100  FORMAT (AS. 151 
IF(jOP  -  1JC«>  (S. 2.0 

K  PAi'SP  h 

7  WPTTr  ri^TP'iT  TAPF  4«ini,TJOP 

101  format  eru  »as//37h  folatton$  rr  mction  fcr  a  towfc  Ronv//> 

TF'sr**sF  Si-'lTfH  ?.)■- Ann,4Pl 

aoo  of  ad  i  ►*Pi*T  tapf  7 .991  .••xsFr.NOOF.  in 
091  format  him 

FAT  Fll  F  10 

OF  AC  I'Pl'T  TAPC  3,097.  (  F^T  (  T  1  .  T  =  1  «  1  7  ) 

097  FORMAT ( 1 7AA ) 

of  a  o  impmt  TAPFi,9o?.  n°f  n  »t«if  n  tiSAdi  .men  ,i*i  ,5» 

092  FORMA  T  (  t  7  !  M 

PFAO  1  f 'Pi  <T  TAPF  3,FMT.(SF(n  *TfOMCn  .1*1.51 
AO l  OFAD  INPUT  TAPF  3.120.T*'fMt*OA.TMWPY,TyWPZtTMFWiWR»Yl»Y2«Y4.Tl  TP*TL 
1C0.TI  Ml  tTfJ»TMtTJ*CMViCtlM«rMP*Tr'1 
120  FORMAT (4F15.Q) 

PFAD  INPUT  TAPF  7,99(7, TNOP,TI  m.TNX,Fy.F7 
990  FORMAT (IA.4F9. 4) 

OFAD  INPUT  TApp  3.999«»'JPA 
999  FOPMA  T  f 1 4 ) 

OFAD  INPUT  TAPr  3.99F *TMuT.T|  T.rjP.TjNP.OT 
99B  FORMAT (5F10. 71 

00  ?P0  I  a  l.NJOF 

OFAD  Input  TAPF  3.120.  ( *'T  ( YK  1  ,KK  si.  4) 

DFAD  INPUT  TAPF  3.170.  (TrMFG(YY) .  YYsl.4| 

OFAD  INPUT  TApr  3. 120. JTTMf (J) .TTMMC7 <j) .TTMMfY(.l) »J*I *5) 

OFAD  JNpi-T  TAPF  3,170.  (TTYOM(.ll  ,TTKCV(.n  .  J«1»?0J 
OF  AC  INPUT  TAPF  3.  120, (TTyTCJI .TTKH(J) .J*1.4) 

OFAD  If.'Pi'T  TAPF  7,120,  (TrOT<J)  ,.|el,4J 
t  S'  ;M  s  C 

NO  ?04  Ktc s  1  ,4 

V  3  -TUCK) 

ryrc A  =  TC*eo.  (YK) 

?K~  =  WTI«1 
tat  =  T  T Y  T  ( <7  ) 

T<M  =  TTYPfYY) 

"0  ?0I  I  s  1,5 
Tm-  =  TT»C(LI 
TM>jC7  s  TTMHC2  fl  ) 

TM«-CY  s  T T'-'HOY  T I  ) 

TS'  V  s  T  stjM  ♦  I 
TY"«  sTTKCMdSUM) 

TK^V  s  T TKCV  (  I  SllM ) 

TMTY  r  TVO  ♦  TvhPY  ♦  TmFw  ♦  T«r/7.  ♦  TVNCY/3. 

TMT7  =  TMP  ♦  TNH07  ♦  T*«Fu>  ♦  TMMT  ♦  TmC/7.  ♦T**wr7/3. 

TLr  T*  s  (  Y1*T"P  ♦  Y2*T,,HP7  «  y4»TMFv»  *TI  TUMWT  *dl  TD/3.  J»(TV£  . 
lT.yiC7J )  / T '■  T 7 
T<.  s  v  R  *  T N 
A  s  TI.PBVH 

Tf-73  (T.lO  ♦Tjwpr  /  e  A  *  Tl  n  1 


COMPUTER  PROGRAM  (Cont’d) 

TIY  «  (TKCM*TH)/A 

TI7  *  (TKCVMTLCM  -TlTPM/A 

TIC  *(  TKCH*TH»«2  ♦  TKCV*(TLCM  -  TLTP)**2  ♦  TKT*fTI  T  -  TLC*>  *T!CD 
l*TO  ♦TK8«<TLCM  -  TLRL)  ♦  T*W) / (A*TLB1 
TTY1  *  (CBH*TJ) /A 

TI7i  a(CBV*(TLCM  -TLCDl  ♦  CT*(TLCM  -  TLT) ) /A 

TIC1  *  (C8H*TJ«*2  ♦  CRV* (TLCM  -  TS.CD>**2  ♦  CT*(TLCm  -  TLT)**2  ♦ 

I  CRR  1/(A*TLB) 

TITY2  *  TMTY/Wfl 

TIIY  *  TKCH/WB 

TITO  *  CT1CCH*TM  -  TKM/A 

TTTY1  *  CBH/WR 

Tirol  *  (CfiH*TJ)/A  ♦  CBP/(A«TlR» 

TIIIZ2  *  TMT7/WR 
Till?  a  TKCV/W8 

TIT  10  a  (TKCV*(TLCM  •  TLT01  ♦  TKB  •  TKT»/A 
TIIIZ1  «  (CBV  ♦CTl/WR 

TIT  101  *  (  CBV* (TLCM  -  TlCM-  CT«(TLT-  TLCM)  ♦  C8R/TlB)/A 

PCAP  a  (  TMX  ♦  TKCV* (TLCM  -  TlTP)*W)/A 

YCAP  *  FY/WB 

ZCAP  «  (FZ  ♦  W*TKCV)/WR 

IF <  MODI)  ?0l*7O0»T'0i 

70!  AA  a  TI IY1*  (  TIT  IZl *TI01  •  T!IT01*TIZl)  •  TnOl*lTTYl*T!llZn 
CYV  *(TIIIZ1*TI01-  TI I IQl *TIZl ) /AA 
CZY  «(-  TlII01*TtY1)/AA 
COY  «(TIY1*T!II7!)/AA 
CY7  *  (-TIIG1*  TIZDMA 
CZ7  a(TIOl*TIIYl  -  T!lOt*TTYl)/AA 
CQ7  *(-TIIYl*TIZll/AA 
CYO  *  IT! I01*TI I IZl ) /AA 
C7"  a{„TIIYl*  TI  I IOH  /AA 
COO  *  <TIIYl*TIIIZ!)/AA 

RYY  a  CYY* (TIIY  -  nVEGA**2«TI IY2 )  -  TIY«CrO 

RYZ  *  CYZ* (TI I IZ  -  0MEGA*«?*T!IIZ2)  ♦  CYQ*T?Z 

RYO  *  CY7»TI!I0  -  CYY«TI!0  ♦  CYO*(TIO  -TIQ2«OMFGA**2) 

RZY  i  CZY* (TI IY  -  TI IY?*0MFGA**2)  -  TIY*CZO 
RZ7  *  CZQ*TI7  ♦  CZ7* { T T I T Z  -  TIT IZ2*0MFGA**2) 

RZO  a  -C7Y*TI 10  ♦  CZZ*TlIIO  ♦  CZCMTIC  -  TIQ2*0MFGA**?> 

ROV  a  CCY»(TTIY  -  TI  TY?*PMFGA**?)  -  COO«TIY 

R07  *  COZMTIIIZ  *TI I I72*QMEGA**2)  ♦  CQC*TIZ 

PQO  *  -C0Y*TII0  ♦  CQ7*rI  T 10  ♦  CCO*(TIQ  -  TI07*0MFGA**?) 

HYY  *  (  RYY**7  ♦  RYZ*R7Y  ♦RYQARCY) /0MEGA**2  ♦  I. 

UZY  «  («ZY*RYY  ♦PZ7*RZY  ♦R70*R0Y1 /0MFGA«*2 
HQv  a  (RQY*RYY  ♦RCZ*RZY  ♦R'.0*R0Y)  /OMEGAS 
HYZ  *  (RYY*RYZ  ♦RZ7*RY7  <-  ROZ*RYO> /0MEGA**2 
HZZ  *  (RZY*RYZ  tRZZ**?  *ROZ*RZO ) /0MEGA**2  ♦  1. 

HQ7  =  (RYZ*ROY  ♦  R77*ROZ  ♦RQQ*PC7 > /0MEGA**2 
HYC  *  (RYY*RYC  ♦  RYZ*R?Q  ♦  ROC*RYQ) /QMEGA#*2 
HZO  a  (RYO»RZY  ♦  P7C*R7Z  ♦  RZG*ROC) /0MFGA**2 
HOC  *  (RYORQY  ♦  ROQ*«?  ♦  RQZ*R7Q) /0MEGA*#2  ♦  1. 

BR  a  HYY»(H7?*HQ0  -MZO*HOZ)  ♦  HYZ*(HZO#HQY  -WZY*HOQ)  ♦  HYQ*(HZY» 

1 HQZ  -  HZZ*HQY) 

TJYY  a (HZZ*HQO  -NZO*H07  ) /RB 
TJZY  «  (HZQ*HQY  *  HZY*HQO)/SB 

TJOY  a  (MZY*H07  -HZ7*H0Y)/RR 

TJYZ  «  (HYQ*HG?  -  HY7*WG0)/RR 

TJ7Z  «  (HYY*HCO  -  HYO*»JQY)/BB 
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COMPUTER  PROGRAM  (Cont'd) 

TJGZ  •  IHY2«HQY  *  MYV«HC7)/B8 
T.JVO  •  (HVZ*H7G  -  HZ7«WVnj/B8 
TJZO  «  IHYO«HZY  *  HYY«HZOJ/BB 
TJCO  *  fHYY«H7Z  -  hV7«h7y1/BB 
E  •  CYY»YCAP  ♦  CYZ*ZC*P  ♦  fYG«Pf*P 
F  a  CZY«YC*P  ♦  CZZ»ZfAo  ♦  CZQaPCAP 
6  a  COY*VCAP  ♦  COZ*ZCAP  ♦  COO«PfAP 
Y1  *  (TJVY*E  ♦  TJYZ*F  ♦TJY0«P1/C«F6A 
Z1  a  (TJZY*E  ♦  TJ7Z«F  ♦TJZG*GJ/CKEGA 
01  *  (TJOY*E  ♦  TJ OZ«F  ♦TJGQ«P  1/OPEGA 
Y2  *  (RYY»Yt  ♦  RYZ*Zl  ♦  PYG*01  1 /O^EGA 
77  *  (R7Y«Y1  ♦  Q?Z«71  ♦  o?C»C!  1 /C*E6A 
07  *(  RCY*Y1  ♦  RC7*71  ♦  PGO«Ctl 7C*EGA 

710  YA*'P  s  5QRTFf  Yt  **?  ♦  Y7»«21 
7A**P  s  S0RTE(71##7  ♦  7?**?> 

PA-P  s  SPPTF(C1**Z  ♦  C7**?t 

tf f scnsf  switch  31  900,901 

9C0  WMfTP  OUTPUT  TAPE  4 , 1 04 ,TI T V7 .Tf JI 7? .T 107 ,TT f Y1 «Tf TG! ,TI 1 171 *TI HC 
U.-nYltTIZl.TI01.TUY.TlT0.T!II7.TinG«TIY«TI7.Tl0.YCAP»ZCAD,PCAP 
1 04  FO"MATMM  «  AW  TIIV2  a. FIS. A,  OH  *11122  s*Fl5.f>,7H  Tt02  a«F1S.P, 

1  AH  rt*VS  a.FtS.ft.  AH  TlTfl  »»F15.°*/  9W  TIIlZl  a,F15.P*  9H  T  T I IQ1 
7a.?15.».7H  TtYl  a,F1S.*».7H  TI71  a.F15.P.?W  TI01  r.FtS.B./  7W  TIlY 
Ss*ri«;.A.7h  TUG  =. FIS. a. AM  TIIT7  s,F!5.°*BH  Til!"  s,Fi5.3,6H  TIY 
4«F}S.8/6l>  TI?  s.F!*;.qfAW  T?0  =«F1S.8,7M  YCAP  *,FlS.B,7h  *CAP  =  ,F1 
5.B.7H  PCAP  s.FIS.8///) 

<s01  wrtTIF  OUTPUT  TAPE  4,151 

IS!  FORMAT  (1  MO,  11 7M  *«HT  I.T  CT  J8  JHP  OVFC.A 

1  ST  *P  MC  :JvC7  mhcy  YCW  KCV  W 

75 

•JRTTP  OUTPUT  TAPE  4,1  S^tT^MT ,TI  T»CT »TjP*TJHPtnMFGA,TYT,TYB,T''C,TMW 
lfZ*TMHCY»T<CW,T£CV*W 
ISO  FORMAT (1  MO, I? F«. ?, ?F1 0.7, FA. ?) 

!F{M"D!)  703,704,703 

704  CA)  S  ACTM(Yt.Y7,TMF7A,n 
YA*'GI  =  THETA 

CAM  ArTM7l»7?.TMFTA«T) 

ZA“GI  s  TWFTA 

CAM  ACTN(C! ,0?,TMFTA»T1 

OA“GL  a  THETA 

Y3  a  -Yt«Cn5F(YAAGl  )  -  Y?»SIf'EfYANG|  ) 

73  s  -71*COSc<7AfcGi )  -  7?*SIHF (7A6.'Gl  ) 

03  a  -ri*COSF  (OAFGl  )  -  p?*MMf  fCANGL) 

GO  TO  705 

703  CAM  ACTf)(V2»Yl»TMFTA,n 
YA*’GL  a  THETA 
CAM  ArTM(72,71,THFTAtn 
7A*  61,  r  THETA 
CAi  L  ACTf  (0?»ni  ,TmFTA.T) 

OA'-GI  =  THETA 

Y3  a  -Y 1  *S I NF  ( vAf'Gl  )  -  V7«C0SF  fyAA'^l  } 

Z3  s-Zt«SlUF(7ANGL)  -  72*CO?F (7A”6l) 

03  S  -01  *Sp'F  (OAA'Gt  J  -  '■‘?«rnSE(0AV6n 

705  IF (  Y3)  604,407,607 

60S  IF  f 73)  604*607, fn7 

60 A  IF (C1)  609*607,607 

601  VPtTF  01  Ti'l'T  TAPE  4*6’<o,  ■n*71iC.> 

410  FORMAT  II"  ,  1  Aw  TRCCO  »’OT  maxI'V*  3*70.3) 

609  WRrTF  OUTPUT  T a Of  4,lS6*Yt  ,Y?,7i ,77,0i,C2,YAMP,?AMpt0A^P»YANGI  *ZAN 
t  GL .C  AROL 
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COMPUTER  PROGRAM  (Cont'd) 


FORMAT CT HO.  SH  Y!  *,F1«:.a«  5H  Y7  s«F15.A*5w  71  «»F15.*«5H  22  *»F15*8 
1,8.  5H  Cl  *.F15.P«  5H  02  *.F15.8//7h  YA«P  «»F15.8«7H  7AMP  «»F15.8» 
?7h  rAmp  *»F15.A//8h  VA»6l  *»F15.8«9H  7Af.GL  »»F15.8«8M  CAijGL  *«F15» 
181 

IF I SFN5F  SWITCH  1)  6?0.*?l 
820  Yin  «  -TTY 
const  s  n 
Y01  r  -TTY1 
01  IDT  s  -TT 101 
01  If)  *  — T ITC 

601  FCPMAT (3F20.8) 

WRITE  OUTPUT  TAPF  6.A0T#TIO2»TICl»TietCCNST.Yni.YID«CnNST*TI71*TI2 

i  .const .pcap. const 

WRITE  OUTPUT  TAPr  A .603, CONST .01 101  tOl  !C «T  1 1 Y?  .Tl  I Yl  «T 1 1 YtCONST *C0 
INST. CONST »CONST.YCAP. CONST 

WRITE  OUTPUT  TAPE  6,603«r0NST»Tl 1 101 ,TI I IC.COMST .CONST *CCNST,TI I !2 
12«TTTl21,TIII7.CONAT«2rAP,fONST 
621  IF  (SENSE  SWITCH  7)602.701 

602  CAl L  SFTHPI  IR <11 «5 , I S ( 1 1 • TO (1) • I SS C 1 » *  TO) 

TINC  *  I NCR 

TINCO  *  1.0/TTNC 

TMAX  s  6.281853/OMFGA 

IJ  *  TNAX/TINCR  ♦  U 

00  995  JJ  a  1.IJ 

TI  *  JJ  -  1 

ANG  «  TI«TINCR*OMFGA 

Y  s  V1#C0SF(ANC.)  ♦  Y?*«;lMF(AN6J 
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Appendix  B 

COMPUTATION  OF  HYDRODYNAMIC  MASS,  DAMPING. 
AND  SPRING  CONSTANTS 


COMPUTATION  OF  HYDRODYNAMIC 
MASS  CONSTANTS 

As  shown  in  Eq,  (11),  the  mass  matrix 
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must  be  evaluated.  If  the  body  has  a  plane  of  symmetry  defined  by 
gravity  and  the  fore-and-aft  direction,  and  if  the  coordinate  axes 
are  along  the  principal  axes  of  inertia,  this  matrix  in  its  complete 
form  would  be 
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where 

mt  is  the  tota!  mass  in  a  given  direction. 


The  total* mass  is  composed  of  the  body  mass  (or  rotational  inertia), 
the  flood-water  mass  (or  rotational  inertia),  and  the  hydrodynamic 
mass  (or  rotational  inertia).  The  cross-coupling  "masses,"  mtya  , 


mfz-  ,  m,-y  and  mtaz  ,  are  hydrodynamic  masses  only, 
because  of  the  above  choice  of  coordinate  axes. 

Hydrodynamic  masses  for  bodies  of  various  shapes  moving  in 
translation  have  been  determined  experimentally  by  Patten.81  The 
effects  of  frequency  and  displacement  amplitude  have  been  investi¬ 
gated  by  Miller.  *  Hydrodynamic  mass  moments  of  inertia  can  be 
computed  for  ellipsoids  in  accordance  with  the  methods  of  Zahm.83 
Information  can  not  be  found  in  the  literature  that  enables  one  to 
compute  the  cross-coupled  hydrodynamic  masses.  Consequently, 
the  mass  matrix  used  in  Eq.  (12)  has  all  off-diagonal  elements 
equal  to  zero. 

COMPUTATION  OF  DAMPING  CONSTANTS 

It  is  extremely  difficult  to  find  data  in  the  literature  that  enable 
one  to  compute  the  damping  constants  for  a  towed  body.  Newman®4 
presents  an  analytical  technique  to  compute  damping  coefficients  of 
ellipsoids.  Patton®5  presents  experimental  data  for  model  towed 
bodies  moving  in  translational  motion.  These  data  can  be  scaled 
to  a  full-size  body  if  scaling  laws  for  size  and  frequency  can  be 
determined.  This  method  was  used  in  the  computation  of  damping 
constants  for  the  AN/SQA-11  VDS  towed  body.86  It  was  assumed 
that  the  damping  varies  with  the  cube  of  the  body  size. 
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Experimentally  measured  rotational  damping  constants  and 
cross-coupled  damping  constants  can  not  be  found  in  the  -literature. 
The  damping  o f  the  AN/SQA- 1 1  towed  body  oscillating  in  the  pitch 
mode  was  approximated  by  computing  the  damping  on  an  equivalent 
ellipsoid  after  the  method  of  Newman, 


COMPUTATION  OF  SPRING  CONSTANTS 

The  spring  constants  due  to  the  body  lift  and  the  tail  lift  have 
been  considered  in  the  text.  A  certain  amount  of  judgment  is  called 
for  at  this  point  because  the  springs  may  be  nonlinear.  (AC(k 
and  AC may  be  nonlinear.)  One  must  estimate  the  pitch 
amplitude  and  substitute  the  best  approximate  linear  values  for  the 
actual  nonlinear  values. 

The  cable  spring  constants  and  Kcv  can  be  computed 

by  considering  two  springs  in  series.  The  first  spring  represents 
the  axial  extension  of  the  cable  under  some  load.  This  spring  con¬ 
stant  is  given  by 


AE 
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where 

A  =  cross-sectional  area  of  the  cable, 
E  =  cable  modulus  of  elasticity,  and 
L  =  length  of  the  cable. 


The  second  spring  in  the  series  combination  is  due  to  the  changes 
in  configuration  of  the  towcable  because  of  changes  in  load  at 
the  towed  body.  First,  the  equilibrium  configuration  of  the  system 
is  computed  by  using  suitable  tables  or,  possibly,  the  computer 
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solution  programmed  by  Cuthill.**7  The  drag  of  the  body  is  then 
increased  by  some  amount  typical  of  the  dynamic  loads  existing 
in  the  cable.  The  new  configuration  is  computed,  and  the  change 
in  trail  distance  is  used  to  compute  a  spring  constant  K  yy  .  The 
change  in  depth  can  be  used  to  compute  a  cross-coupled  spring 
constant  Kyz  .  The  same  process  is  undertaken  to  compute  the 
spring  constants  Kzz  and  K  zy  by  increasing  the  weight  of 
the  body. 

Equivalent  spring  constants  and  KCT  are  computed 

as  follows: 
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ch 


KyyKlS»»9c 

Kyy  +  Kx  sin  d0 


and 

K  =  K«K1co,<&b 
CV  =  K^  +  Kjcos^’ 

where 

<f>0  is  the  angle  (from  the  vertical)  of  the  towline  at  the  towed  body. 

The  spring  constants  Kc|,  and  Kcv  are  not  constant 
because  the  springs  that  they  represent  are  nonlinear.  One  should 
compute  a  number  of  deflections  for  different  increments  of  drag 
and  weight  to  approximate  the  nonlinear  springs  in  the  yy-  and 
zz-directions  with  linear  spring  constants  Kclj  and  Kcv  . 


b7H.  Cuthill,  A  FORTRAN  Program  for  the  Calculation  of  the  Equilibrium  Configuration  of 
a  Flexible  Cable  in  a  Uniform  Stream,  David  Taylor  Model  Basm  Report  No.  182,  May  1963 
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